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The production of a DD meson-pair in antiproton-proton (pp) annihilation close to the production 
threshold is investigated, with special emphasis on the role played by the '!/>(3770) resonance. The 
study is performed in a meson-baryon model where the elementary charm production process is 
described by baryon exchange. Effects of the interactions in the initial and final states are taken 
into account rigorously, where the latter involves also those due to the '!^(3770). The predictions for 
the DD production cross section are in the range of 30 - 250 nb, the contribution from the '!/i(3770) 
resonance itself amounts to roughly 20 - 80 nb. 
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I. INTRODUCTION 


In a recent publication [l| we have presented predic¬ 
tions for the charm-production reaction pp DD close 
to the threshold based on a model where the elemen¬ 
tary charm production process is described by baryon 
exchange and also in the constituent quark model. The 
cross section was found to be in the order of 10 - 100 nb 
and it turned out to be comparable to those raedicted by 
other model calculations in the literature [2 t3- 

The results in Ref. [ij suggested that the reaction 
pp —> DD takes place predominantly in the s wave, 
at least for excess energies below 100 MeV. However, 
there is a well-established p-wave resonance, the '0(3770) 
= 1 ) which is seen as a pronounced structure 

in e+e” —>■ DD ii , for example, and which is located 
at only around 35 MeV above the DD threshold. The 
resonance decays almost exclusively (i.e. to 93^9 %) into 
DD 0- This resonance was not included in our pre¬ 
vious study [Ij. Given its apparent prominence in the 
DD channel, the impact of the '0(3770) on the pp —>■ DD 
cross section clearly should be explored. In particular, 
should it turn out that its contribution is rather large, 
then it would be certainly interesting to examine the en¬ 
ergy range in question in pertinent experiments, which 
could be performed by the PANDA Collaboration [n^ 
in at the future FAIR facility in Darmstadt. 

In the present study we consider the effect of the 
i/'(3770) on the pp —)> DD cross section. The work com¬ 
plements our results presented in Ref. [l| and builds on 
the Julich meson-baryon model for the reaction pp —>■ 
KK [I^ where the extension of the model from the 
strangeness to the charm sector follows a strategy sim¬ 
ilar to other studies by us on the DN and DN inter¬ 
actions [IMl, and on the reaction pp A^ A+ [T^ , 
namely by assuming as a working hypothesis SU(4) sym¬ 
metry constraints. 


II. THE MODEL 

The framework in which we treat the charm- 
production reaction pp DD is described in detail in 
Ref. [l|. Here we summarize only the principal features. 
The reaction amplitude for pp —>■ DD is obtained within 
the distorted-wave Born approximation. Effects of the 
initial as well as of the final-state interactions, which 
play an important role for energies near the produc¬ 
tion threshold MM. are taken into account rigorously. 
The employed NN and DD amplitudes are solutions of 
Lippmann-Schwinger type scattering equations based on 
corresponding {NN and DD) interaction potentials. 

The microscopic charm-production process itself is de¬ 
scribed by baryon exchange (Ac, Sc), in close analogy 
to an investigation of the strangeness moduction reac¬ 
tion pp KK by the Julich group [IJ. Specifically, 
the transition potential is derived from the correspond¬ 
ing transition in the KK case under the assumption of 
SU(4) symmetry, see Ref. [I| for details. 

Because of the known sensitivity of the results for 
the cross sections on the initial pp interaction we ex¬ 
amined its effect by considering several variants of the 
NN potential. Details of those potentials can be found 
in Refs. UM- Here we just want to mention that they 
differ primarily in the elastic part where we consider vari¬ 
ations from keeping only the longest ranged contribution 
(one-pion exchange) to taking a full G-parity transformed 
NN interaction as done in [l^. All those models repro¬ 
duce the total pp cross section in the relevant energy 
range and, in general, describe also data on integrated 
elastic and charge-exchange cross sections and even pp 
differential cross sections, cf. Refs. (ll. [l8l|. 

The interaction in the DD system is constructed along 
the lines of the Julich meson exchange model for the tttt 
interaction whose evaluation has been discussed in de¬ 
tail in Refs. [U, The present interaction is based 
on the version described in the latter reference. The po¬ 
tentials for TTTT TTTT, TTiT ^ KK and KK ^ KK are 
generated from the diagrams shown in Fig. [TJ The addi¬ 
tional diagrams that arise for the DD potential and for 
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FIG. 1. Diagrams included in the Jiilich tttt — KK potential 
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FIG. 2. Additional diagrams that arise when the DD channel 
is included. 


the transitions from tttt and/or KK to DD are shown in 
Fig. 3 In this extension we were guided by SU(4) symme¬ 
try Thus, we included t-channel exchanges of those 
vector mesons which are from the same SU(4) multiplet 
as those included in the original Jiilich model and, more¬ 
over, we assumed that all coupling constants at the addi¬ 
tional three-meson vertices are given by SU(4) relations. 
Indeed, in Ref. [I| an even more extended model was con¬ 
sidered which included also the coupling to channels in¬ 
volving the charmed strange meson IIs(1969). It turned 
out that the DD production cross sections based on the 
DD interactions with or without coupling to D'^D~ are 
almost identical and, therefore, in the present work we 
show only the results for the latter case. 

The scattering amplitudes are obtained by solving a 
coupled channel scattering equation for these potentials 
which is formally given by 

i 

with i,/, Z = TTTT, TTT], KK, DD. 


In an approach like ours the ^/>(3770) (in the following 
usually called ip to simplify the notation) has to be in¬ 
cluded as bare resonance. It aquires its physical mass and 
its width when the corresponding potential is iterated in 
the scattering equation ©• We include the ip only in the 
direct Uli potential, i.e. . The corresponding 

diagram is depicted on the lower right side of Fig. [2] The 
potential can be written in the form 


yDD.DD ^ ^DD 


1 


E — mo 




( 2 ) 


with a bare DDip vertex function and a bare 

'0(3770) mass toq. Explicit expressions for the vertex 
function of two pseudoscalar mesons coupled to a vector 
meson in the s-channel and for the resulting potential can 
be found in the Appendix of Ref. . The bare mass toq 
and the bare coupling constant in ^o^ are adjusted in 
such a way that the resulting T-matrix, , has a 

pole at the physical values of the 0(3770) resonance. The 
main results shown in the present study are based on a 
DD model that produces a pole at E — (3773 — i 13.6) 
MeV, i.e. at the value specified as “our fit” by the PDG 
0 , but we consider also the value obtained in Ref. @ 
for the mass which is about 6 MeV larger. In the ac¬ 
tual calculation the values for the mass and width of the 
0(3770) were determined by evaluating the speed plot for 
rpDD.DD £gj. sijnpiicity reasons. In addition, we use an 
isospin averaged mass for the D (D), namely 1866.9 MeV. 

Cross sections for DD scattering in the isospin 1 = 0 
and 1 = 1 states can be found in Fig. |31 For the 
1 = 0 case we show the result from the s wave sepa¬ 
rately (dashed curve) so that one can see the impact due 
to the 0(3770) resonance. Note that there is actually a 
noticeable non-resonant contribution in the p wave which 
is clearly visible for energies away from the 0(3770), i.e. 
around 3.9 GeV and above. The cross section in the 
J = 1 state is very small. It is only in the order of 1 mb 
and, therefore, hardly visible in Fig. |31 In this context 
let us mention that there are other model results on the 
DD interaction in the literature , achieved likewise in 
a meson-exchange approach. 

For completeness reasons we present here also results 
for the reaction ^ DD, cf. Fig. |4l Results cor¬ 

responding to a 0 mass of 3773 or 3779 MeV, respec¬ 
tively, are indicated by the thick and thin solid lines. 
The pertinent calculation was performed in the Migdal- 
Watson approximation [^, i.e. by assuming that 
rpDD,e+e ^ yDD,DD^ cross section for 

e“'"e“ DD is given by 


^e+e-^DD 


NqnD 


(3) 


Here qjyp) is the DD center-of-mass momentum and N 
an arbitrary normalization factor that has to be adjusted 
to the data. The shown curve is based on the DD 1 = 0 
amplitude in the p wave. In this case a factor has 
to be divided out because in e+e" —?> DD the incoming 
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FIG. 3. Total cross section for DD scattering as a function 
of the center-of-mass energy for isospin 1 = 0 (solid line) and 
7 = 1 (dash-dotted line). The s-wave contribution to the 
7 = 0 cross section is indicated by the dashed line. 



FIG. 4. Integrated cross sections for e^e“ —>■ DD. Our 
results based on the Migdal-Watson approximation are shown 
by the thick solid line (m.^ = 3773 MeV) and the thin solid 
line (m^ = 3779 MeV), respectively. Data are taken from 

Ref. 1. 


777? state is not on-shell and, therefore, does not feel the 
angular-momentum threshold factor |^ . 

Our result is in remarkably good agreement with the 
measured line shape. In particular, the strong fall-off 
after the peak value is very well described. As mentioned, 
we adjusted our resonance parameters to the “our fit” 
value of the PDG. No fine tuning to the actual data was 
done because we are primarily interested in studying the 
effect of the tj;{3770) onpp —>■ DD observables and not on 
pinning down its resonance parameters. As already said, 
there is an uncertainty in the order of 6 MeV with regard 
to the actual resonance mass 0,13. For more detailed 
analyses of the e+e“ —>■ DD data, see Refs. 0,0,[2^[23- 

m. 


III. RESULTS FOR THE REACTION pp DD 

Results for the reaction pp —>■ DD are presented in 
Figs. [S] and IHl We performed again calculations for all 
pp ISI considered in our previous work [l|. However, 
we refrain from showing bands here because we wanted 
to highlight the impact of the ip{3770) and this infor¬ 
mation would have been otherwise hidden in the band. 
Nonetheless, the variations due to the ISI can be still 
read off because we present curves corresponding to the 
largest and smallest predictions for the pp 77°7)° and 
pp —>■ 77+77” cross sections. 

The solid and dash-dotted lines in Figs. [5] and |6] are 
results with inclusion of the ■0(3770) resonance while the 
dashed and dotted curves are based on the DD FSI with¬ 
out the 0, i.e. correspond to the case considered already 
in Ref. [J. Please recall that in our calculation the 0 is 
only included in the 7777 interaction and, therefore, its 
effects come exclusively from the corresponding FSI. In 
principle, there could be also a coupling of the (bare) 0 
to the NN system which would then contribute to the 
pp —>■ DD transition potential. However, we assume here 
that this possible coupling is negligibly small - because 
it is suppressed by the OZI rule. After all, it requires the 
annihilation of 3 up- and/or down quark pairs and the 
creation of a cc quark pair; see discussion in Sec. El 

It is obvious that the 0(3770) resonance leads to a pro¬ 
nounced enhancement in the 7777 production cross sec¬ 
tion. Especially, for the 77° 77° final state the cross section 
predicted around the resonance peak is now practically 
twice as large. Quantitatively, our calculation suggests 
that the cross section due to the 0(3770) resonance could 
be the order of 20 to 80 nb. The effect is less dramatic in 
the 77+77” case because here the s-wave contributions are 
fairly large. Still for both reactions there should be good 
chances that one can see a clear resonance signal over the 
background (that comes mostly from the s wave) in an 
actual experiment. 

A comparison of the two results with and of the two 
without 0(3770) resonance (solid versus dashed line or 
dash-dotted versus dotted line, respectively) in Fig. [5] al¬ 
lows one to see the uncertainties in the predictions due 
to the employed pp interactions. Those are substantial 
but not really dramatic as already discussed in Ref. [I|. 
There is a somewhat larger uncertainty with regard to the 
contribution of the 0(3770) resonance which, however, is 
not too suprising. Conservation of total angular momen¬ 
tum and parity implies that the 7777 s wave is produced 
from the NN ^Pq partial wave while the DD p wave is 
produced from the NN ^Si-^Di partial wave. Since the 
partial wave is much more sensitive to short-range 
physics than higher partial waves there is also a stronger 
variation in the corresponding amplitudes for the various 
NN potentials considered as ISI. 

The thin solid lines indicate what happens if we employ 
a 7777 FSI that is fitted to the higher value given for the 
0(3770) resonance, i.e. to 3779 MeV 0,[i3- In this case 
the contribution of the 0 to the pp —>■ DD cross section 
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FIG. 5. Total reaction cross sections for pp —>■ DD as a function of piab- The (red) solid and dashed curves and the (black) 
dash-dotted and dotted curves are results for different pp initial-state interactions, see text. Solid and dash-dotted curves show 
the full results with inclusion of the 'i/>(3770) in the FSI, while dashed and dotted curves are without i/)(3770). The (magenta) 
thin solid curves indicate results based on a DD FSI fitted to the higher tp mass (3779 MeV) obtained in Ref. [^. 




FIG. 6. Differential cross sections for pp —>■ DD at piab = 6.578 GeV/c (excess energy e = 40 MeV). Same description of curves 
as in Fig. [5] 


is smaller by about 15-20 % and, of course, the peak is 
shifted to the higher energy. 


IV. ESTIMATES OF THE CROSS SECTION 

FOR pp —>■ DD 


Differential cross sections for pp —>■ and pp —>■ 

D^D~ at Piab = 6.578 GeV/c are shown in Fig. [SI This 
momentum corresponds to a total energy of 3773.8 MeV 
so that the results display the situation practically at 
the resonance peak. As mentioned before [T| without 
the ^/'(3770) the DD pair is produced predominantly in 
the s wave, cf. the dashed/dotted curves. Of course, 
this changes drastically once the resonance is included 
(solid/dash-dotted curves). Still, there is a strong in¬ 
terference between the s- and p waves, especially in the 
D^D~ case, so that the actual angular distribution does 
not resemble the one one expects from pure p-wave scat¬ 
tering. 


It is possible to provide a rough estimation of the 
pp —>■ DD cross section around the '0(3770) peak using 
experimental information by exploiting the fact that at 
energies very close to the resonance the reaction ampli¬ 
tude can be well approximated by 


rj~\l 


E — -I- ir,/,/2 


(4) 


an expression which is indeed exact at the '0(3770) pole. 
Here, 7 ® is the dressed vertex function with i,j = DD, 
pp, e+e”, and and F,^ the physical mass and width 
of the 0(3770), respectively. 

The cross section resulting from this amplitude can be 
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cast into the form 

(2J+1) 47r r,/2 rj/2 

(25f + 1)(25'’ + 1)^{E- + ry4 ’ 

(5) 


i.e. into the standard Breit-Wigner formula. Here J is 
the total angular momentum, S'® and are the spins of 
the particles in the incoming channel, qi is the center-of- 
mass momentum in the incoming channel, and Ti and Fj 
are the partial widths for the decay of the '!/'(3770) into 
the channels i and j. 

There are cross section data on the reactions e+e” —> 
DD and e+e“ —>■ pp [3l| in the vicinity of the 
'!/'(3770) resonance. Together with unitarity constraints 
for DD —7> DD they allow one to determine all the T^’s 
needed for estimating the pp DD cross section. With 
regard to the DD —>■ DD cross section, we assume for 
simplicity that Tjjq ss T^ - which is anyway consis¬ 
tent with the PDG listing 0 - so that the cross sec¬ 
tion is given by the unitarity bound. At the energy 
corresponding to the ^’(3770) resonance this amounts to 
o'dd^dd ~ 200 mb. Assuming that the 1 = 1 p-wave 
amplitude is small, which is indeed the case in our model 
(cf. Fig. [3]) yields ~ 50 mb. The mea¬ 

sured e+e“ —5> D^D~ cross section at the '!/’(3770) reso¬ 
nance peaks at about 3 nb @ while the analysis of the 
e+e“ —)• pp cross section performed in Ref. suggests 
that the contribution due to the ■!/'(3770) resonance could 
amount to either 0.059to.o2o 2.57to;i3 pb. Combining 

those results and using Eq. to determine the various 
Ti’s we deduce for pp D~^D~ a cross section of either 
3 ^ 18 nb or around 350 nb. There is a large uncertainty 
for the lower value that follows from the BESIII analysis 
but, still, one could argue that it is roughly in line with 
our theory prediction. 

As alternative let us consider also an estimate that fol¬ 
lows from a QCD-based perturbative approach to char- 
monium decays into baryon-antibaryon pairs [s^, [s^. 
That approach relies on a factorized expression for the 

^ pp decay amplitude into a cc ^ pp annihilation am¬ 
plitude and the ip wave function at the origin - see, e.g., 
Eqs. (4.68)-(4.71) of Ref. [S^- The model dependence 
regarding the cc —>■ pp annihilation amplitude, which in¬ 
volves at least three gluons in the creation of up and down 
quark-antiquark pairs, can be eliminated considering the 
ratio of T^^pp to Tj/^^pp. If, in addition, one assumes 
that the strong as well as the electronic decays of ip{3770) 
proceed predominantly through the 2 component of 


its wave function (recall that the existence of an impor¬ 
tant 2 component comes from the large decay width 
of ip{3770) into e+e“), one can express F^^pp in terms 
of measured quantities: 


^ _( I- \ " Fy,^e+,- 

1 - • 

Using PDG values [10|| for the electronic decay widths 
and for the width of J/dt —>■ pp, one obtains from Eq. ([S]) 
for the cross section <Jpp^D+D- ~ 0.2 nb, due to the 
■0(3770), which is significantly lower than the value ex¬ 
tracted from the BESIII result but also much smaller 
than our model predictions. 


V. SUMMARY 

We have presented a study of the reaction pp DD 
close to the production threshold with special emphasis 
on the role played by the 0(3770) resonance. The work 
is an extension of a recent calculation by us [I| which is 
performed within a meson-baryon model where the ele¬ 
mentary charm production process is described by baryon 
exchange, and where effects of the interactions in the ini¬ 
tial {pp) and final {DD) states are taken into account 
rigorously. The 0(3770) resonance, which is included in 
the DD ESI in the present calculation, produces a size¬ 
able enhancement in the pp DD cross section around 
the resonance energy. Indeed, our predictions for the to¬ 
tal DD production cross section in the considered near¬ 
threshold region are in the order of 30 - 250 nb, where the 
contribution from the 0(3770) resonance itself amounts 
to roughly 20 - 80 nb. Given the magnitude and the 
shape of the cross section due to the 0(3770) there should 
be good chances to measure the pertinent contribution 
in dedicated experiments which could be performed at 
FAIR. 
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